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Summary 


This report provides an analysis of experimental laser Doppler 
velocimetry (LDV) data uncertainties that propagate from measurements in 
the tunnel coordinate system to results in the model system. Calculations of 
uncertainties as functions of the variables that comprise the final result 
requires assessment of the contribution each variable makes. Such an 
analysis enables and necessitates the experimentalists to identify and 
address the contributing error sources in the experimental measurement 
system. This provides an opportunity to improve the quality of data derived 
from experimental systems. This is especially important in experiments 
where small changes in test conditions are expected to produce small, 
detectable changes in results. In addition, the need for high-quality 
experimental data for CFD method validation demands a thorough 
assessment of experimental uncertainty. 

Transforming from one Cartesian coordinate system to another by 
three sequential rotations, equations were developed to transform the 
variables initially obtained in the original coordinates into variables in the 
final coordinate system. Based on the transformation equations, propagation 
equations for errors in the experimentally-derived flow quantities were 
derived for a model at angle of attack. Experimental uncertainties were 
propagated from the tunnel coordinate system into the model system. 
Comparisons between results for the two systems revealed a variety of 
increases and decreases in bias and precision errors. 
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Introduction 


Changes in experimental test conditions, with expected changes in 
test results, require high repeatability of experimental data and estimates of 
experimental uncertainties. For example, the validation of computational 
fluid dynamics (CFD) methods requires experimental data of high quality 
with the best estimates of experimental uncertainties. Laser Doppler 
velocimetry (LDV) is frequently used in wind and water tunnels to make 
measurements in off-body flow fields. This technique provides accurate 
measurements of flow mean and fluctuating velocities from which various 
turbulence relationships, e.g., Reynolds normal and shear stresses, can be 
calculated. 

Many publications of experimental LDV results include an analysis of 
expected errors or uncertainty in some form. Several documents are 
dedicated specifically to error or uncertainty analysis. References 1-3 are 
examples of such analyses. Experimental results are often transformed into 
coordinate systems other than those in which they were measured. For 
example, if one or more of the axes in a multi-axis LDV system are not 
congruent with the tunnel axes, the measured quantities must be resolved 
into the tunnel system. In addition, if the model being tested is at an angle 
relative to the tunnel coordinate system, the measured quantities must also 
be transformed into the model coordinate system. 

The propagation of experimental uncertainties through the various 
coordinate transformations is required to provide accurate estimates of 
uncertainties in the final, transformed experimental results. Systematic 
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procedures for propagating uncertainty have been described by several 
authors. One of the most current references in this area is by Coleman and 
Steele (reference 4). 

This report addresses the propagation of experimental uncertainties 
in LDV flow field measurements from the tunnel coordinate system into the 
model coordinate system. The purpose is to provide a systematic way to 
transform contributing errors of various calculated flow quantities into final 
results. 
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Symbols and Abbreviations 


a, A 
ANSI 
ASME 

b, B 
B 

c, C 

d, D 
e 

LDV 

N 

P 

r 

RSS 

S 

t95 

U 

u,v,w 
U, V, w 
u’.v'.w' 

Vu ’ 2 , Vw i2 

u'u', v'v', w'w' 
uV, v'w', w'u' 


transformation matrix elements 
American National Standards Institute 
American Society of Mechanical Engineers 
transformation matrix elements 
bias limit 

transformation matrix elements 
transformation matrix elements 
unit vector 

laser doppler velocimetry 
number of samples 
precision limit, P = tgsS 
result in an uncertainty analysis 
root-sum-square 
precision index 

statistical t-distribution value giving 95% 
confidence level 
uncertainty or velocity, in/ sec 
total instantaneous velocities, in/ sec 
mean velocity components, in/sec 
instantaneous velocity fluctuations, in /sec 
rms velocity fluctuations, in /sec 
Reynolds normal stresses, in 2 /sec 2 
Reynolds shear stresses, in 2 /sec 2 
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X 


x, y, z 
a 

P 

A 

5iJ 

4 > 

P 


variable in uncertainty analysis 
cartesian coordinates, in 
angle of attack, deg 
sideslip angle, deg 

implies "uncertainty" when preceeding 
quantities in uncertainty analysis 
Kronecker delta 
roll angle, deg 

correlation coefficient between 
variables in uncertainty analysis 


subscripts: 

B u 

i or j 

ij 

j 

LDV 

p ij 

r 

rms 

RSS 

t 

x 

x, y, z 

oo 

1 . 2 , 3 


correlation coefficient for bias errors, eq (2) 
elemental component 
matrix element indices 
jet 

LDV coordinate system 
correlation coefficient for precision errors, 
equation (2) 
result 

root mean square 
root-sum-square 
tunnel 
variable 

unit vector axes 

free stream condition 

intermediate axis systems 
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superscripts: 


-> 

A 


fluctuating quantity 
mean quantity 
vector 
unit vector 


Analysis of Experimental Uncertainties 

The concepts of bias and precision errors are the fixed and random 
errors, respectively, that occur in experimentation. The bias errors are 
fixed, systematic or constant errors that induce an offset from the true value 
of the quantity being measured. The precision errors are random variation 
or repeatability errors of the quantity being measured. The total uncertainty 
is then formed by a combination of a bias limit and a precision error 
estimate. The bias limit, B, is an estimated limit "of a confidence interval on 
the true value of the bias" (reference 4). The combination technique used in 
the present study is the root-sum-square (RSS) method, where the total 
uncertainty is given by, 


U rss =(B 2 +P 2 x ) ,/2 , 

and where B is the bias limit estimate and Px is the precision limit. The 
precision limit is the product of the precision index (precision error 
estimate) and tgs, the value of t from the statistical t distribution that gives a 
95% confidence level to the estimate of the precision (random) error. For 
experimental sample sizes greater than 30, the standard in reference 5 
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recommends using a value for tgs of 2 for 95% confidence estimates. The 
equation for the total uncertainty becomes, 

Urss = [ B 2 + (2S) 2 ] 1/2 , (1) 

where S is the precision error estimate or precision index. 

Propagation of Uncertainty 

Coleman and Steele (reference 4) describe the propagation of 
uncertainty in a general uncertainty analysis. If a result, r, is a function of 
many variables (xi, X2, X3, . . . , xn), the uncertainty, U r , in a calculated result 
is 
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The Ux’s are the uncertainties in the variables, xj. In a detailed uncertainty 
analysis, the bias and precision limit estimates must be propagated 
separately. Then the estimates are combined into the total uncertainty 
(using the RSS method in this case). If the bias limits (or precision limits) of 
the different variables are not independent of each other, there are cross 
terms in the expressions for the bias and precision limits as follows. 

j_ 

2 

, and 


B r = •s 


1=1 


dr 

dx, 


B, 


| ^ dr dr 
+ £ dx, dXj 
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P r = <! 


1=1 



+ f.^^P„P,,Pp,(l-8,) 


pdx l 0Xj 


(2) 


i 

} 

I 

! 

\ 


where the p's are correlation coefficients associated with each type of error 
and 5y is one if i = j and zero if i * j. If the error terms are independent (i.e., 
the errors in the measurements of the variables are independent), the cross 
terms dissapear (p = 0) in equations (2). 


The bias limit for r is 


B r = 


,3x, x ' 


_§Lr 

3x„ B x3 


_*L b 

3 x n x » 


J 


and the precision limit for r is 


P r = 


— p 
3x, x 


V f 
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* P ' 


3Xn 


f dr _ 
+ ■•■+ - — P x 

V^Xn 
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* 



where the B x 's and the P x 's are the bias and precision limits for the 
measured variables, xj. The equation for the propagated precision error 
estimate (precision index) is 


S r = 


3r 

3x, 


iL s ) 

’-J 


(JLs ) 

dx„ - I 


V 


J 


where the Sx's are the estimated precision errors in the measured results. 
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In a data acquisition program, the results for the measurement 
instrument may be modified by at least one coordinate transformation to 
yield final results. For example, there may be a transformation from the LDV 
coordinate system to the tunnel coordinate system, if one or more of the 
LDV axes are not aligned with the tunnel axes. Then, if the model is at an 
angle or angles relative to the tunnel, there is a transformation to the model 
coordinate system. Therefore, the uncertainty analysis is required to account 
for this transformation. The propagation of the uncertainties into the 
calculated results, including coordinate transformation, was carried out in a 
manner similar to that used by Neuhart, et al (reference 6), where the 
transformation from the LDV coordinate system to the tunnel system was 
developed. 

The description of the development of the coordinate transformation 
equations, from the tunnel to the model coordinate system, are described in 
the following section. Next, the method for estimating uncertainties in the 
final, transformed results will be shown. Finally, an example calculation is 
given to demonstrate the process. Uncertainties in the calculated quantities 
from reference 6, in the tunnel coordinate system for a forebody /strake 
model at angle of attack, are listed in table 1 . The results of propagation of 
the uncertainties into the model coordinate system are listed in table 2. The 
tabulated results show the propagated precision error estimates (precision 
indicies) instead of the precision limit. The total uncertainty, however, does 
contain the precision limit, P = tgsS. All quantities calculated were 
normalized by appropriate terms, as shown in the tables. Final estimated 
uncertainties in table 2 are presented with two significant figure accuracy 
since more significant figures would be inappropriate for such estimates. 
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Coordinate Transformation 


The derivation of the transformation equations follows the 
development of Neuhart, et al., (reference 6), which was based on the work 
of Morrison, et al, in reference 7, In figure 1, the x-y-z axes shown form the 
orthogonal, model coordinate system. The total velocity vector defined in 
this system would be 

U = ue x + ve y + we z , 

where the e x , e y , e z vector quantities represent unit vectors in the x. y, and z 

directions, respectively. This represents the velocity defined in the model 
coordinate system. The velocity defined in the tunnel coordinate system (xi, 
yi, zi axes in figure 1) would be represented generally as 

U. = u t e x _ +v t e yi +w t e Zi , 

where the e X| , e X| , e Xi quantities are the unit vectors in the tunnel system 

coordinate directions. In general, the two coordinate systems have the same 
origin, but the axes are not coincident. By successive rotations (similar to 
Euler rotations), the transformation equations between the tunnel and the 
model at sideslip, angle of attack, and/or roll can be obtained. 

Beginning with the xi, yi, z i axes in figure 1, a rotation is made about 
the zi axis (sideslip, (3). As shown in figure 2, the relationship between the 
respective unit vectors can be represented by 
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or in matrix form, 


e x , =e X 2 cosp-e y2 sinp 
e yi = e Xj sin p + e y2 cosp. 


e 

z i 

e z 2 




A’ 


cosp 

-sinp 

O' 

^x 2 

A 

= 

sinp 

cosp 

0 

A 

«... 


0 

0 

1 

A._ 


A subsequent rotation (figure 3) about the y 2 axis (angle of attack, a) yields, 




e v 

= e K cosa + e 7 sina 

* 


x 2 

*3 z 3 



A 

= 



^Z 2 

= -e_ sina + e z cos a 

x 3 z 3 


or, 


A 


cos a 

0 

sina 

V 

A 

=r 

0 

1 

0 

A 

A,. 


-sina 

0 

cos a 

_^ Z 3 _ 


A final rotation about the X 3 axis (roll, <j), figure 4) results in 


e ya = e y cos<()-e z sin(J), 
e Z3 = e y sin<|) + e z cos<J> 


A* 


'1 

0 

0 

A 

A 

= 

0 

COS<j) 

-sin<|) 

e y 

A,. 


0 

sin<)) 

COS()) 

A. 
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The relationship between the unit vectors in the tunnel coordinate system 
(xi.yi.zi) and those in the model system (x,y,z) is, by substitution for 
intermediate unit vectors, the product of the three matricies. 


A 


A' 


'cosp 

-sinp 

O' 

" cos a 

0 

sin a" 

i 

0 

0 

A 

A 

- 

A, 

= 

sinp 

cosP 

0 

0 

1 

0 

0 

COS(|> 

-sin<() 

A 

_K_ 


A_ 


0 

0 

1 

-sin a 

0 

cos a 

0 

sin<J> 

COS<J> 

A. 


Therefore, the relationship between the velocity vectors is 


u, 


'cosp 

-sinp 

O' 

’ cos a 

0 

sin a 

"1 

0 

0 

"u" 

V, 

= 

sinp 

cosp 

0 

0 

1 

0 

0 

cos<|> 

-sin<() 

V 

_ w t. 


0 

0 

1 

-sin a 

0 

cos a 

0 

sin<|> 

cos<t> 

w 


or in shorthand form, 


A ' 

v t 

=N 

[ B »I C »] 

u 

V 

-N 

’u 

V 

. w t. 



w 


w 


Since the desired result will be in the model coordinate system, the 
inverse process will be used. This process is defined by 


“u " 


u t' 


’V 

V 

=[ c d l [B,r l [A 1J ]-' 

v t 

■hr 

v t 

w 


. w «. 


. w t. 


Utilizing the rule that the product of a matrix and it's inverse equals the 
identity matrix, the terms in each rotation matrix were determined. The 
first equation, for sideslip, is shown below. 
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cosp 

-sinp 

O' 

" a n 

a !2 

a i3 


"1 

0 

o' 

sinp 

cosp 

0 

a 21 

a 22 

a 32 

= 

0 

1 

0 

0 

0 

1 

„ a 31 

a 32 

a 33_ 


0 

0 

1 


Multiplying the terms to form nine equations, the ay terms were 
determined. The solution matrix for sideslip, ay, is 

’ cosP sinp 0" 

-sinp cos (3 0 . 

0 0 1 


The above inverse matrix, ay, was verified by multiplying it by the original. 
Ay, matrix to yield the identity matrix. 

The transformation matrix for angle of attack was derived using the 
following equation. 


" cosa 

0 

sina 

'b„ 

^12 

^13 


1 

0 

o' 

0 

1 

0 

^21 

^*22 

^32 

= 

0 

1 

0 

-sina 

0 

cosa 

_^31 

^32 

^33^ 


0 

0 

1 


The solution matrix for angle of attack, by, is 

cos a 0 - sin a" 

0 10 . 

sina 0 cosa 

Finally, the transformation matrix for model roll was obtained using 
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1 

0 

0 

C ll 

C 12 

C 13 


1 

0 

o' 

0 

cos<}> 

-sin<]) 

C 21 

C 22 

C 32 

= 

0 

1 

0 

0 

sin<(> 

cos<|> 

_ C 31 

C 32 



0 

0 

1 


The transformation matrix for roll, cy, is 

10 0 
0 cos<f> sin<|> . 
0 - sin <J> cos 0 


The resulting full transformation from the tunnel to body coordinate system 
is given by equation (3) as 


u " 


'i 

0 

0 ' 

cos a 

0 

-sina' 

’ cos 3 

sinP 

O' 

"U|‘ 

V 

= 

0 

cos<{> 

sin<|> 

0 

1 

0 

-sinP 

cosP 

0 

v t 

W 


0 

-sin<}> 

cos<}> 

sin a 

0 

cos a 

0 

0 

1 

. w t. 


The three separate matricies can be combined to form a single matrix for 
the three body rotations. 


"u " 


cos a cos P 

cos a sin P 

-sina 


V 

w 

— 

sin a cos P sin <|> - sin p cos <J> 
sin a cos P cos <)> + sin p sin <|> 

sin a sin p sin <(> + cos P cos <> 
sin a sin P cos <> - cos P sin <{> 

cos a sin 4> 
cosacos<|> 

v t 

_ w t 


or in short form, 


u 

V 

II 

jET 

i 

i 

< c 

i 

=N 

u t ' 

V t 

w 


L w tJ 


. w t. 


( 4 ) 
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Since the total instantaneous velocities are the sum of a mean and an 
instantaneous fluctuating velocity, 

u = u + u’ 

V = V + v' , 
w = w + w' 

then the following transformations exist. 


u 


’ n t" 

V 

w 

=N 

V, 

_w t _ 


(5) 


u‘" 


'u t ' 

v l 

w’ 


V t ' 

. w t’_ 


From the definition of the Reynolds normal and shear stresses, 


and 


N 

X u i u i 

u'u' =— uu 

N . (6) 

N 

X U i V i 

u'v' = — uv 

N . (7) 


where Uj and vj are toted instantaneous velocity values. Substituting the 
matrix equations (4) and (5) for the ui and u values into equation (6) yields 
the following transformation equation for the u'u' Reynolds normal stress. 
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1=1 1. 

u'u' = — - 


^ll^ll U t U t "*■ ^ll^l2 U l V t ■*" ^11^13 U t W t 

+d I2 d n v t u, +d 12 d 12 v t v t +d, 2 d 13 v t w t 
+d 13 d,,w t u t + d I3 d 12 w t v t +d l3 d 13 w t w t 
N 


1 


^d n d, , u t u t + d, ,d 12 u t v, + d, ,d I3 u t w t ^ 

+d 12 d u ^ + d 12 d 12 v^v^ + d I2 d l3 v;w; 


+d 13 d, ,w t u t + d 13 d 12 w t v t + d I3 d 13 w t w t 


Substituting for the ui, vi, u, and v quantities in equation (7) from equations 
(4) and (5) gives the transformation equation for the uV Reynolds shear 
stress. 


N f dnd 2 ,u,u t + d M d 22 u t v t +d M d 23 u t w t 


i-1 


u 




+d| 2 d 21 v t u t + d 12 d 22 v t v t + di 2 d 23 v t w t 
+d 13 d 21 w t u t +d 13 d 22 w t v t + d I3 d 23 w t w t 


l 


_N 

^d„d 21 u t u t +d n d 22 u t v t + d„d 23 u t w t A 

+d, 2 d 2 , vj’uT + d 12 d 22 + d 12 d 23 v^ 


+d 13 d 2 , w t u t + d, 3 d 22 w t v t + d l3 d 23 w t w t 


Similar types of derivations were done for the other two normal 
stresses and two shear stresses, and are not shown here. Using the rule that 

the summation of a sum of terms equals the sum of the summations of the 
individual terms (^(a + b + c) = £a + £b + £c) and rearranging terms 

yields transformation equations that are a function of normal and shear 
stresses in the tunnel coordinate system. The result is a set of 
transformation equations that are expressed as the following matrix 
equation. 
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u'u’ 


di jdj j 

di 2 di2 

d^dia 

2d u d 12 


2d n d 13 

2di 2 d 13 


V 

v T v’ 


^21^21 

d 22 d 22 

d 23 d 2 3 

2d 2 id 22 


2d 21 d 23 

2d 22 d 23 

v t 

v l’ 

w' w 1 


^31^31 

d3 2 d 32 

d 3 3d 3 3 

2d 3 id 32 


2d 31 d 33 

2d 32 d 33 

W, 

W,' 

u 1 v' 


di id 2 i 

d 12^22 

di 3 d 23 

dnd 22 + d 12 d 21 

d 

nd 23 + d 13 d 2 | 

^12^23 + di 3 d 22 

u t 

’ v t' 

w’ u’ 


daidi i 

d 3 2di 2 

d 33 d i3 

d 3 idi2 + d 32 d n 

d 

3id 13 + d 33 di j 

d 32 d 13 + d 33 d 12 

W t 

'U t ' 

v' w' 


_d 2 id 3 i 

d 22 d 32 

d 23 d33 

d 2 id 32 + d 22 d 3J 

d 2 id 33 + d 23 d 31 

d 22 d 33 + d 23 d 32 

_v 

w t 'J 


In an uncertainty analysis, the contribution of the bias and precision 
limits of each variable determines the total uncertainty of the result. Clearly, 
each resultant Reynolds stress error term can be can be a function of up to 
nine variables. These are the six tunnel system stresses plus the three 
rotation angles in the dy matricies. For simplicity, we will consider the most 
common case where a # 0 and (5 = <|> = 0. The resulting transformation matrix 
for the mean velocities is 



cos a 0 - sin a 
0 1 0 
sin a 0 cos a 


and for the Reynolds normal and shear stresses, 


u'u’ 


cos 2 a 

0 

sin 2 a 

0 

-2 cos a sin a 

0 

u t 'u t ’ 

v'v' 


0 

1 

0 

0 

0 

0 

W 

w'w' 


sin 2 a 

0 

cos 2 a 

0 

2cosasina 

0 

w t ’ w t ’ 

u'v’ 


0 

0 

0 

cos a 

0 

-sin a 

u ’v ‘ 

Ll t V t 

w'u* 


cos a sin a 

0 

-cos a sin a 

0 

cos 2 a - sin 2 a 

0 

w t ’u t ' 

_v'w’ 


0 

0 

0 

sin a 

0 

cos a 

.v t 'w t ' 
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The resulting equations for the transformed mean and instantaneous 
fluctuating velocities are 


and 


u = u t cosa- w t sina 

v = vj" 

w = u t sina + w t cosa 

u ' = u t 'cosa- w t 'sina 
v’= v t ' 

w'= u t 'sina + w t 'cosa 


( 8 ) 


The transformation equations for the Reynolds stresses are 

u'u' = u t ' u t ' cos 2 a + w t ' w t ' sin 2 a - 2w, ’ u t ’ cos a sin a 
v'v' = v t 'v t ' 

w'w' = u t ' u t ' sin 2 a + w t ' w t ' cos 2 a + 2w, ' u, 'cos a sin a 


( 9 ) 


uY = u t ' v t 'cosa - v t 'w t 'sina 

w'u' = u t ' u t ' cos a sin a - w, ' w t ' cos a sin a + w t ' u t '(cos 2 a - sin 2 a) ►. 
v'w' = u t 'v t ' sin a + v t ' w t 1 cos a 


( 10 ) 


The transformed rms fluctuating velocities Uj-ms. v rms. and w rms are 
calculated as the square root of the transformed u'u', v'v', and w'w', 
respectively, after transformation. This is done since, for example, 



(ID 


19 



Also, therefore, 


V — a/v' v'. 

v rms v v » 

w = Vw' w' 

V¥ rms ’ ¥¥ 


The transformation equations for the rms velocity fluctuations involve the 
square root of transformation equations (9) for u'u', v'v’, and w'w'. 

These equations contain Reynolds normal and shear stresses in the tunnel 

coordinate system (see equations (9)). However, according to equation (11), 
Ut'u,' = u tn J, and also v t 'v t ' = v t ^ 2 and w t 'w t ' = w,^ 2 . Therefore, using 

equations (9) and (11), 


=• 

Vu'u' =( 

u. 2 cos 2 a + w t 2 sin 2 a-2w.'u,' 

ran l rm* 1 1 

cosasina) 1/2 

V ™ S =‘ 

\Zv'v’ = y 

'v.'v,' 


W rms = 

V w'w' = 

fu. 2 sin 2 a + w t 2 cos 2 a + 2w,'u, 

\ rma 1 rms l 1 

"T . \ 1 / 2 

t cos asm a j 


( 12 ) 


Propagation of Uncertainty into the Model Coordinate System 

Refering to the transformation equations (8) for the mean velocities, 
the equations for the propagation of the bias limits and precision indices 
(errors) can be derived using equations (2). Since 


and 


v = v 


t * 



and, for example, 
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I2 


B- = 

V 


' B-' 


Sv, 


the propagated bias limits and precision errors for the v mean velocities are 
equal to those for the measured quantities. 


That is. 


B- = B— , and S- 



For the u and w velocity components, the bias limit equations become (from 
equations (2)) 


+ f4^B 

U w t w, j 

du du „ D 

+ au^3w^ p “'"’' * s * 


" ^3U, 


' 2 


— \2 


=B- 
3a a 


and 


\2 


B - 2 = 

w 




3w r f 3 w 
3u t Ut J [dw t W *J 


3w 3w _ D 
+ 3n;3^; pu ' w * ^ ^ 


2 X-.— 


^ 3a 


■)’ 


The bias limit in a was ±0.226° (±0.004 rad). The assumption is made that 
the errors in the measurements of the variables ut, wt, and a are 
independent. The cross terms drop out (p UjWt = p Ut „ = p w<a = 0), resulting in 

(from equations (8)) 

B- 2 = (cosa B-) 2 +(-sina B— ) 2 +(-u['sina-w^cosa) 2 B 5 2 and, 

B- 2 = (sin a B-) 2 +(cosa B— ) 2 +(u^cosa-w t sina) 2 B- 2 . 
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Similarly, for the precision errors, 

S- 2 = (cosa S— j +|-sinaS— j .and 
S- 2 = (sin a S— j + (cosa S— j . 


The uncertainty in the angle a is limited to the bias limit only. 


As before, the total uncertainty is then calculated as 



The transformed rms fluctuating velocities were calculated as the 
square root of the transformed Reynolds normal stresses (equations (12)). 
Therefore, similar to the mean velocities, the bias and precision errors in 
v rms quantities were 



The u and w rms fluctuating velocities were given by the transformation 
equations (12). Cross terms are provided for the 

u tnn> and w t 'u t ', and w tmj and w t 'u t ' pairs since they are related in the 
w t ’u t ' covariances. Assuming errors in measurements of 
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u tnns , w lnna , Eind a , were independent (eliminating those cross terms), the 
equations for the bias limits were derived from 
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The assumption is made that the elements in the cross terms are perfectly 
correlated, giving correlation coeffcients of one. This assumption, yielding 
conservative uncertainty estimates, will be adopted in all further 
developments. 


Inserting the appropriate terms in the above equations and 
performing some algebraic manipulation yields 


TV 
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B„ 
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The precision errors were defined by 
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The resulting equations for the precision errors in the Urm S and w rms 
quantities are 


( u t cos 2 a „ 
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The bias limits and precision errors for the Reynolds normal stresses 
were calculated using the transformation equations (9). Again, since 

v'v' = V l , V t ' 


the bias limits and precision error estimates for v’v’ were derived as 


B— = B~^— and S— = S-^. 


v'v' V, V 


V V V, V, 


Assuming that the errors in the measurements of u t 'u t \ w t 'w t ', and a 
are independent for this quantity (eliminating those cross terms), the 
equations for the bias limits in u'u' can be derived from 
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Based on the transformation equation (9) for the u'u' Reynolds normal 
stress, the bias limit is 
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B— 2 =(cos 2 aB— — :) +(sin z aB — - — + (-2cosasina B ) 

UU \ u,u,7 \ w,w,7 \ W t 'u ,7 
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The bias limit in w'w' can be formulated from 
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Using the transformation equation for w'w' from equation (9), the bias limit 
is expressed as 

B^-; 2 = (sin 2 a B-r— + (cos 2 aB— — +(2cosasina B — - — 

ww \ u t u, / \ w,-w,7 \ w t 'u t / 

+(2 cos a sin a(u t ' u, ' - w t ' w t ') + 2 w t ' u t '(l - 2 sin 2 a)j 2 B a 2 
+2 cos a sin a ( sin 2 a B— ;B — — + cos 2 a B— T — ,B^^ ) 

V U, U, W, U, W( W( w t 'u,7 


The precision indicies had a similar form, without the a terms, given as, 


S— r-; 2 = (cos 2 ot S — \ — ;) +(sin 2 aS— 7 — + (-2cosasina S — —) 

u u \ u t 'u,7 \ w,'w, / \ w t 'u,7 
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\ u t u, w, u, W, w, w, u, ) 


and, 
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S^? 2 = (sin 2 a + (cos 2 a S^.f + (2cosasina S^f 

+2 cos a sin a(sin 2 a + cos 2 a S^— -S— - ) 


The uncertainty analysis for the Reynolds shear stresses was based on 
the transformation equations (10). In this case, the LDV velocity 
measurements were required to be made from the same seeding particle for 
correlation of velocity fluctuations. Therefore, the u’v’, and v'w' stresses 
each had a third, cross term in their expressions for bias limits and 
precision errors. The w'u' stresses had two cross terms. The equations for 
determination of the bias limits for the three shear stresses can be derived 
from 
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The equations for the bias limits are as follows. 
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The precision error equations were determined from the following 
relationships. 
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The resulting equations for the precision error estmates are, 
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Example calculation. The calculation of the bias limit in the 


v'w' 


U 


nondimensional shear stress is presented in this section. The equation for 
estimating this error is listed below. 


B-^-t 2 = (sina B— tt) + (cosaB— — =) 

V w \ u t v t / \ v t w t / 

+ (u t ' v t 'cosa - v t 'w t 'sina) 2 B a 2 
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U t V, V, w t 


All velocity quantities used in this equation were nondimensionalized by Uoo 2 , 
yielding the bias limit in the nondimensional shear stress. The angle a is 
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25°. The bias limit in this angle, B a , is ±0.226° (±0.004 rad). The bias limits 

in Ut and Vt were obtained from table lc and are ±0.0000029 and 

Uj 2 Uj 


±0.0000012, respectively. Average measured values of 


U.' V f 

U „ 2 


and 


v t' w , 

XJj 


were 


used and they were 0.0029 and 0.0016, respectively. (If this analysis was 
being done on a point-by-point basis, the values at the particular point in the 
tunnel would be used.) The four terms, in order of their appearance in the 
bias error equation, are 1.50 X 10' 12 , 1.18 X 10 -12 , 6.10 X 10 _1 *, and 1.33 X 
10 ~ 12 . The total bias limit is ±0.0000081. Clearly, the third term dominates 
this error, which is the error in the angle of attack measurement. This 
result is typical of most systems, and better accuracy in angle of attack 
measurement will have the greatest effect in reducing this bias error. 


This concludes the discussion of the propagation of uncertainty for the 
calculated quantities. The tabulated results for all of the transformed, final 
quantities in the model coordinate system are given in tables 2a, 2b, and 2c. 
All transformed quantities were normalized by the appropriate terms, as 
shown in the tables. In all cases where cross terms existed in the error 
equations, the sign of the correlation coefficient was unknown. It was 
therefore assumed to be the sign that would maximize the error estimate. 

Analysis of Uncertainty Results 


The significance of the final uncertainty results can be examined by 
comparing them to values in the tunnel coordinate system. The uncertainty 
results relative to the results in the tunnel coordinate system showed 
changes in all quantities associated with the u and w velocity components. 
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Since the effect of angle of attack, a, was the only effect examined, and the u 
and w vectors only are affected by a rotations, all v-only related quantities 
remained unchanged. All quantities associated with the w-related variables 
showed the largest differences from the values in the tunnel system. 


In general, bias limits and precision errors in the six Reynolds 
stresses increased, except for the errors in 


v'v' 


U 


, which did not change. 


There was a slight decrease in uncertainty in the u rms fluctuating velocity 
due to decreases in bias limit and precision error. The uncertainty in the 
Wrms fluctuating velocity increased primarily due to an increase in bias limit. 
Bias limit increases offset precision error decreases in the mean u velocity. 
Both errors increased for the mean w velocity. 


To address the effect of model sideslip and/or roll on uncertainties in 
the mean velocities, the full transformation in equation (4) would be used to 
formulate the transformation equations. The number of variables in the bias 
error propagation equations would increase by the number of additional 
angular rotations applied to the model. To evaluate the effect on rms 
fluctuating velocities and the Reynolds stresses, sequential application of the 
angular rotations would reduce complexity. Although not proven here, it is 
assumed that this would involve developing the propagation equations for 
each angle separately, setting the other angles equal to zero. The equations 
for each angular rotation would then be applied separately to each 
succeeding result. 
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Conclusions 

An analysis of the propagation of estimated experimental LDV data 
uncertainties from the tunnel coordinate system to the model system has 
been presented. Transforming from one Cartesian coordinate system to 
another by three sequential rotations, equations were developed to 
transform the variables initially obtained in the original coordinates into 
variables in the final coordinate system. Based on the transformation 
equations, propagation equations for errors in the experimentally-derived 
flow quantities were derived for a model at angle of attack. Experimental 
uncertainties were then propagated from the tunnel coordinate system into 
the model system. 

Comparisons between results for the two systems revealed a variety of 
increases and decreases in bias and precision estimates. Quantities 
associated with the w-related variables primarily increased and showed the 
largest differences from the values in the tunnel system. In some cases, u- 
related quantities decreased after transformation, or did not change 
significantly. 

Calculations of uncertainties as functions of the variables that comprise 
the experimental results allows assessment of the contribution each variable 
makes. Such an analysis enables and necessitates the experimentalists to 
identify and address the contributing error sources in the experimental 
measurement system. This provides an opportunity to improve the quality of 
data derived from experimental systems. This is especially important in 
experiments where small changes in test conditions are expected to 
produce small, detectable changes in results. In addition, the need for high- 
quality experimental data for CFD method validation demands a thorough 
assessment of experimental uncertainty. 
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Table la 


Uncertainty Estimates for Results in Tunnel Coordinate System- 
Mean and Fluctuating (rms) Velocities 




(Reference 3 
(Uoo=3 

- Strake Test) 
in/sec) 




u t /U_ 

^/ u „ 

w, /U„ 

V< 2 /u. 

V v t' 2 /U. 

^w,' 2 /U, 

Strake Test 
Bias 

+.0101/ 

+. 0254 / 

+. 00822 / 

±.00195 

±.0119 

±.00441 


-.0112 

-.0281 

-.00957 




Precision 

±.00168 

±.00888 

±.00110 

±.00119 

±.0060 

±.000775 

Total 

+.011/ 

+. 031 / 

+. 0085 / 

±.0031 

±.017 

±.0047 

Uncert. 

-.012 

-.033 

-.0098 
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Table lb 


Uncertainty Estimates for Results in Tunnel Coordinate System 

Reynolds Normal Stresses 
(Reference 3 - Strake Test) 

(Uoo=3 in/sec) 



u t 'u t '/Ul 

v t 'v t ’/U! 

w.’w.’/U; 

Strake Test 

Bias 

0 

±.0000513 

0 

Precision 

±.000056 

±.00149 

±.000024 

Total 

±.000056 

±.0030 

±.000024 


uncert. 



Table lc 


Uncertainty Estimates for Results in Tunnel Coordinate System- 

Reynolds Shear Stresses 
(Reference 3 - Strake Test) 



(Uoo= 

3 in /sec) 



u t 'v t 7u! 

. v t ’w t 7U* 

w t 'u t ’/U 

Strake Test 

Bias 

±.0000029 

±.0000012 

0 

Precis. 

±.00012 

±.000066 

±.000020 

Total 

±.00024 

±.00013 

±.000020 

uncert. 
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Table 2a 


Uncertainty Estimates for Results in Model Coordinate System- 
Mean and Fluctuating (rms) Velocities 




(Reference 3 
(Uoo=3 

- Strake Test) 
in/sec) 

[== 

[== 

Strake Test 

u/U„ 

v/U M 

w/U M 

Vu /U_ 

Vv ’ 2 / u . 

Vw ' 2 / u . 

Bias 

+.0110/ 

-.0120 

+.0254/ 

-.0281 

+.00909/ 

-.0103 

±.00181 

±.0119 

±.00546 

Precision 

±.00159 

±.00888 

±.00122 

±.00117 

±.0060 

±.000771 

Total 

Uncert. 

+.011/ 

-.012 

+.031/ 

-.033 

+.0094/ 

-.011 

±.0030 

±.017 

±.0057 
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Table 2b 


Uncertainty Estimates for Results in Model Coordinate System- 

Reynolds Normal Stresses 
(Reference 3 - Strake Test) 

(U co =3 in/sec) 

tTTT/uE ^v/uE v^/uE 


Strake Test 

Bias ±.0000110 +.0000513 ±.0000110 
Precision ±.000059 ±.00149 ±.0000322 

Total ±.00012 ±.0030 ±.000065 

uncert. 
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Table 2c 


Uncertainty Estimates For Results in Model Coordinate System- 

Reynolds Shear Stresses 
(Reference 3 - Strake Test) 

(Uoo=3 in/sec) 


Strake Test 

uv/ui 

Vw’/Uj, 

w’u'/U* 

Bias 

±.0000111 

±0000081 

±.0000009 

Precis, 

±.000125 

±.0000958 

±.0000339 

Total 

uncert. 

±.00025 

±.00019 

±.000068 
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